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Abstract 

Based on two-photon exchange interaction between n coupled optical cavities each of 
them containing a single three level atom, the n-qubit and n-photonic state transfer is 
investigated. In fact, following the approach of Ref.[l], we consider n coupled cavities 
instead of two cavities and generalize the discussions about quantum state transfer, pho- 
ton transition between cavities and entanglement generations between n atoms. More 
clearly, by employing the consistency of number of photons (the symmetry of Hamil- 
tonian), the hamiltonian of the system is reduced from 3 n dimensional space into 2n 
dimensional one. Moreover, by introducing suitable basis for the atom-cavity state space 
based on Fourier transform, the reduced Hamiltonian is block-diagonalized, with 2 di- 
mensional blocks. Then, the initial state of the system is evolved under the corresponding 
Hamiltonian and the suitable times T at which the initially unentangled atoms, become 
maximally entangled, are determined in terms of the hopping strength £ between cavities. 

Keywords: coupled cavities, two-photon exchange, hopping strength, three 
level atoms, generation of entanglement, excitation and photon transfer, 
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1 Introduction 

The quantum communication between several parts of a physical unit, is a crucial ingredient 
for many quantum information processing protocols [2]. Schemes for the transfer of quantum 
information and the generation and distribution of entanglement have been designed and im- 
plemented, in the past years, in a number of physical systems (see for example [3]- [12]). Atoms 
and ions are particularly considered as tools for storing quantum information in their internal 
states. Naturally, photons represent the best qubit carrier for fast and reliable communication 
over long distances [13, 14]. Recently, using photons in order to achieve efficient quantum trans- 
mission between spatially distant atoms has considered in several works [1, 15, 16, 17, 18, 19]. 
The basic idea, is to utilize strong coupling between optical cavities and the atoms. On the 
other hand, due to the ability of quantum entanglement as a resource for several quantum 
information processing tasks such as quantum communication, and certain quantum crypto- 
graphic protocols, the creation of quantum entanglement naturally arises as goals in nowadays 
quantum control experiments in studying the nonclassical phenomena in quantum physics. 

One of the known models in quantum optics describing the atom-field interaction is the 
Jaynes-Cummings Hamiltonian [20, 21]. In the study of three-level atoms, M. Alexanian and S. 
Bose [17] introduced a unitary transformation, whereby the three- level atom was reduced to a 
corresponding two-level atom of the Jaynes-Cummings type with two-photon instead of single- 
photon transitions. In Refs. [1, 18], entanglement properties of two and three atom-cavity 
systems in which the cavities are coupled via two-photon exchange interaction, was analyzed in 
detail. Such results could set the pathway towards massively correlated multiphoton nonlinear 
quantum optical systems [22, 23], which are rapidly developing modern subjects nowadays. 
The motivation of interest to such systems is their promise in quantum switching, quantum 
communication and computation and quantum phase transition applications. 

In this paper, following the approach of Refs. [1, 17, 18], and introducing some suitable basis 
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in which the Hamiltonian of the system can be block-diagonalized, we generalize the discussions 
of Ref.fl] and Ref.[18] for two and three atom-cavity systems, to a system consisting of n 
coupled atom-cavity subsystems. More clearly, we consider a system of n spatially separated 
optical cavities, each containing a single three level atom, which are coupled to each other 
with two-photon exchange interaction. Our objective is to examine state transfer (atomic 
state exchange or photon transition) within photon and atom subsystems and to consider 
possible generation of the particle entanglement between the subsystems. 

The organization of the paper is as follows. In section 2, the model describing a system 
of n identical atom-cavity subsystems is introduced. The main results of the paper such 
as block-diagonalization of the Hamiltonian of the system, solving the Shrodinger equation 
for time dependent probability amplitudes of the state of the system, and discussions about 
state transfer (atomic excitation or photon transitions) and entanglement generation between 
atoms or photons, are given in this section. Sections 3 and 4 are respectively concerned with 
the special cases of two and three identical coupled cavities. Paper is ended with a brief 
conclusion. 



2 The Model: n coupled cavities via two-photon ex- 
change interaction 

We will consider n identical cavities each containing one three-level atom, where the cavities 
are coupled via two photon hopping between them. In fact, we consider that the cavities are 
located at the nodes of the complete graph K n with n nodes and each cavity interacts with all 
of the other cavities via two-photon exchange. 

Let us first introduce the two-photon Hamiltonian obtained via an exact unitary transfor- 
mation introduced in Ref. [17]: 
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where, the operator 

ATM = a\a % + a« " ^ + 1 

is a constant of motion for the i-th atom-cavity subsystem, i.e, we have [H^\N^] = for 
each i = 1,2, ...,n. The operators a, and a| are photon operators of the i-th cavity, and 
a ab = f° r * = 1,2, ... ,n denote the atomic transition operators for the i-th cavity 

referring to either the ground (g) or excited (e) state. Now, the Hamiltonian for the n cavities 
is given by: 

n n 

H = J2(H^-Ht) + H E {a?a) + afa1), (2-1) 

i=l i,j=l;i<j 

where, 

HP=hu>(N®-l) + (E g + E e )/2, 

with E g , E e being the energies of the ground and exited states, respectively. The last term 
in the Hamiltonian (2-1) is the two-photon exchange interaction between the cavities, char- 
acterized by the hopping rate £. The parameters E ,/i,r) and A are the free energies of the 
subsystems written in the notation of Ref. [17] and we do not need their clear definitions in 
the present paper. All of these parameters depend on the photon number in the corresponding 
cavities and so, on the cavity-mode intensity through the eigenvalues of the operator N®. 

The operator N = Yh=\ commutes with the Hamiltonian (2-1) and so we can reduce 
the Hamiltonian to the subspace spanned with the eigenstates of N and consider the time 
evolution of the states in this subspace. For a given eigenspace of N with eigenvalue N, the 
maximum possible number of photons in a cavity is N when the corresponding atom is in the 
ground state, which occurs when there are no photons present in the other cavities and the 
atoms are also in the ground state. Then, the total number of photons in the system will be N. 
The constant number of total photons determines the subspace or the manifold in which the 
states evolve in time (the initial state of the system determines the constant number N). We 
will consider the manifold with N = 2. In this case, each single atom-cavity system can take 
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one of the three possible states \g,0), \g,2) or |e,0), and so, the total possible states that the 
system of n-cavities can take, are 3™ states. Due to the consistency of total N = 2, the only 
possible states which we can have, are 2n states instead of 3 n ones. In fact, these 2n states 
are eigenstates of N with eigenvalue 2, and the 3 n -dimensional Hamiltonian H is reduced 
to 2n- dimensional one in the bases which span the eigenspace of N with the corresponding 
eigenvalue 2. The bases states that span this subspace or manifold, are given by: 

\d) = | ff ,O)...|0,O)|£^|0,O)...|0,O), 

i—th 

\a l ) = \g,0)...\g,0)M\g,0)...\g,0), 

i—th 

for i — 0, 1, . . . , n — 1. Indeed, these bases span the eigenspace of N with eigenvalue 2, i.e., we 
have N(a\ci) + /?|aj)) = 2(a|cj) + /3|aj)). Therefore, the general time dependent state of the 
n-cavity system is given by 



n— 1 



\m) = J2(^(t)\c t ) + Mt)\a t )). 



(2-2) 



i=0 



Then, one can easily show that, by considering the order of bases as |c ), |a ) 
the Hamiltonian H takes the following direct product form 



I c n-1/) 



H = I n 



tan 9n 



\ 



tan 9 tan 2 9 



+ 2£(J n - I n ) 







(2-3) 



where, I n is n x n identity matrix and J n is the all one matrix of order n. The quantity 
tan^o is given by tan# = with r = where g\ and g 2 are the atom-photon coupling 
constants in the three-level atom. In writing the above equation, the dimensionless time 
[(Eq — Eq) cos 2 9 ]t/h — > t and dimensionless hopping constant H^/[(Eq — Eq) cos 2 9 ] — > £ 
have introduced (see Ref.[l, 18] for the cases n = 2 and n = 3 coupled cavities), where Eq 
and Eq are eigenvalues associated with eigenvectors \4>o)^ = sin^ |e, 0) + cos9 \g,2) and 
\ipo = cos^ole, 0) + sin^olfi',2) of respectively. 
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It is well known that the matrix J n has eigenvalues 0, and n (due to the fact that J 2 
nJ n ), and is diagonalized by discrete Fourier transform F defined as F H := -j^u kl for k,l 



0, 1, . . . , n — 1, where u = exp 



new Fourier transformed bases {|cj)', |aj)'}™ =0 as: 



is the n-th root of unity. Therefore, by introducing the 

n n-1 



n-1 



l«>':=-7=E<A*>> 



1 



71-1 



(2-4) 



i=0 



and considering the ordering { |co)', |ao)'; • • • ; |c n -i)', |a n _i)'}, the Hamiltonian (2-3) takes the 
following block diagonalized form: 



H = L 



tan # 

2 



tan # tan 9 



+ 2£ dia^(n-l, -!,...,-!) 







(2-5) 



where, diag(n — 1, — 1, . . . , — 1) is the n x n diagonal matrix with diagonal entries as n — 1 
and —1 respectively. Now, by using the Schrodinger equation of motion ih-j^\ip) = H\tp), the 
equations of motion are given by: 



iC' Q = [1 + 2£{n - 1)]C + tanMi; 
L4q = tan6» Co + tan 2 9 A' , 



and 



c (0 



^ = (1-20(7; + tan MI; 

iA{ = tan e Q C[ + tan 2 M{ , (2-6) 

for / = 1,2, — 1. The equations (2-6) can be exactly solved for any value of tan# . 

Here we take the ratio of atomic couplings in the three level atoms as r = l/y/2 so that we 
have tan# = 1- Substituting tan# = 1 in (2-6) and solving the corresponding differential 
equations, one can obtain 

e -i[l+£(n-l)]t 



:{[ V /T+ ?(n - 1)2 costal + e(n - iy-ii(n-l) sint^/l + e 2 (n-l)2]C;(0)- 
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e -i[l+£(ra-l)]t 



;sin^l + e(n- 1)^(0)}, 



(i) = {[^1 + £ 2 (n - l) 2 costal + e 2 (n - l) 2 +^(n-l) sin^l + ^ 2 (n - 1) 2 K (0)- 

^/l + £ 2 (n- l) 2 



isin^l + e(n-l) 2 C^(0)} (2-7) 
where, for I — 1, 2, . . . , n — 1 we obtain 

C M = ^^{[^^ 2 ^st^lTe + i ^mt^lTe]C'M -*sintyTT?4(0)}, 

= ^^{[\/^costV / r+?-^sintV / l + e]^(0) -iamty/T+eCim- (2-8) 

By using (2-4), one can obtain the time dependence of the coefficients Cj(i) and A^t) of the 
state of the system in (2-2) via the inverse Fourier transform as, 

1 71—1 

Q(0 = ^=£^'(0, 

V n i=o 
i n—1 

^(*) = -^E w_,< ^(0- (2-9) 
It should be pointed out that, one can evaluate the probabilities associated with the 
state of the system as a superposition of atomic states |a») , and that of photonic states 
denoted by P a (i) and P c (t), respectively. For instance, considering the initial state 
1^(0)) = ^(\g,2)\g,0)...\g,0) + \g,0)\g,2)\g,0)...\g,0) + . . . + \g,0)...\g,0)\g,2)), with initial 
conditions Ai(0) = and C/(0) = ^= for all 1 — 0,1, ...,n — 1, with the aid of Eqs. (2-2) and 
(2-8), we obtain 

n— 1 n—1 i 

Pc(t) = E IQ(0I 2 = E |C;W! 2 = , 1+ 2(n _ 1)2] {e 2 (n - l) 2 + cosHy/l + 



=0 1=0 



n — 1 



-{£ 2 + cos 2 tvTT?}, 



n(l + e) 

n—1 n—1 

p a (t) = E I4(*)l 2 = E I4(*)l 2 = i - Pc(t). (2-10) 

2=0 2=0 

where, in the second equality in P c {t) and that of P a (t), we have used the fact that the Fourier 
transform is unitary and so dose not change the norm of vectors. The above result indicates 
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that, in the limit of large £ — > oo, we have P c (t) ~ 1 for every time t, i.e., for large enough £, 
all of the atoms will be at their ground state \g) at every time t. 

2.1 Large and small hopping strengths 

One should notice that for large values of the hopping strength, i.e., £ ^>, the evaluated 
coefficients C'^t) and A[(t) in (2-7) and (2-8) take the form 

C' (t) ~ e-^^^-^lHe^^-'^C^O) - iSi ^ ( "~ 1) ^ , o(0)}, 
C7/(f) ~ e-^{e^(0) - ^4(0)}, 

4(0 * e-^ie-^A'^O) - ^C;(0)}; Z = 1, 2, . . . , n - 1. (2-11) 
Neglecting also the second terms in the above approximations, we get 

C' (t) « e-^-^C^O), 

C;(t)^e 2i ^(0), / = l,...,n-l, 
w^(0), / = 0,l,...,n-l, 

and so by using (2-9), we obtain 

p 2i£t n—1 n—l 7i—l p 2i£i n— 1 

C7,(f) « — {e" 2 ** £ C fc (0) + u;^]C fe (0)} = — ^(e" 2 ^ - 1 + n5 kJl )C k (Q), 

n k=0 k=0 i=l n fc=0 

« A,(0); for I = 0,1,..., n- 1 (2-12) 

where, in the first relation we have used the fact that for the n-th root of unity u>, we have 
Yh=q u>( k ~V l = n8 k ,i and so Y^=i ^ k ~ l " >l — n$k,i — 1- The above results, are in correspondence 
with those of Refs. [1, 18] for the special cases n = 2 and n — 3. Moreover, the relations (2-12) 
indicate that in the limit of large hopping strength, the state associated with the initially 
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unentangled atoms, i.e., the initial state with Ai(0) = 0, for / = 0, l,...,n — 1, remains 
effectively unentangled forever. 

In the limit of small hopping £ <C, the equations (2-7) and (2-8) lead to the following 
coefficients C[(t) and A'^t) 

C' (t) ~ e -*[i+«(»-i)]*{ [cost - if(n - 1) smt]C' {0) - % sinL4' (0)}, 

A'o(t) ~e- i[1+ « (n - 1)]< {[cost + ^(n- l)sint]^(0) - isintC£(0)}, 

C[{t) ~ e-^-^{[cost + if (n - 1) sint]C;(0) - isinfcA{(0)}, 

~ e- i(1_€) *{[cos* - if (ra - 1) sin^O) - i sintC;(0)}; I = 1, 2, . . . , n - 1. 

(2-13) 

Now, by neglecting the terms proportional to £, the above approximations read as 

C£(t) « e-^+^-WicostC'oiO) -ismtA' Q (0)}, 

A' Q (t) w e-^+^-^^cos^^O) -isin^(O)}, 
C/(f) « e- i(1 -^{costC;(0) - isinL^(O)}, 

« e- i(1 -^{costA;(0) -isintc;(0)}; Z = 1, 2, . . . , n - 1. (2-14) 
Then, by using (2-9), one can obtain for I — 0, 1, . . . , n — 1 

-i(l-£)t n-1 

C,(i) ~ ^(e-^ ni - 1 + ra<S fc) ,)(costC fc (0) - isinL4 fe (0)), 

n fc=0 

~ 6 °* E^"^ " 1 + n5 k>l )(co8tA k (0) - isin«7 fc (0)), (2-15) 

n fc=0 

It could be noted that for times such that £nt <C 1, also for times such that £t = ^ L 
with k G Z, the above result leads to Q(t) = e-^-^cosfC^O) - ismtA^O)) and A,(t) ^ 
e-^-^coatAiiO) -ismtQ(0)), so that we have \Q{t)\ 2 + \Mt)\ 2 = \Q(0)\ 2 + |^(0)| 2 and 
so, there is no exchange between the cavities. On the other hand, for the times such that 
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£t = ( - 2; ^ 1 - >7r , with I £ Z, the exchange between the cavities (excitation or photon trans- 
fer) can be achieved. For instance, in the case of two cavities n — 2, for the initial state 
|^(0)) = \a ) = \e, 0)\g,0) with initial conditions A o (0) = 1 and C o (0) = A^O) = d(0) = 0, 
by using (2-15), we obtain at times t ~ C (t) = A (t) = 0, d(t) = -ie'^ 1 '^ sint and 

A^t) = e-^-^cost, so that we have \if>(t)) = e^ 1 '®* (cos t\g, 0)|e, 0) - ismt\g,0)\g,2}). 

The results of this section can be used in order to discuss about qubit state transfer, photon 
transition and entanglement generation between the atoms. In order to clarify that, how one 
can discuss these arguments, we will consider the special cases of two and three identical 
cavities in the next sections in details. 

3 Two coupled cavities: the case n = 2 

For two cavities (n — 2), by using the relations (2-7)- (2-9), one can calculate 
C ' + C[ e~ u 

V2 ~ vr+e 



C (t) = -^T^ = ^==={[^l + ecos^tcost^l + e-^m^tsmt^l + e]C (0)- 



i [ y/lT? sin £t cos ty/l + £ 2 +£ cos £t sin tyfl + £ 2 ] d (0) -% sintyl + £ 2 (cos &A (0) -% sin %bA x (0) ) } , 



Ci(t) = < ^f L = + ^ sinet costv /i + ^ + ^cos^sini^/l + £ 2 ]C (0) + 

[^1 + £ 2 cos f i costal + £ 2 -£ sin f i sint^/l + £ 2 ] C x (0) -i sintyl + £ 2 (-i sin £L4 (0)+cos £L4i(0))}, 



4,(t) = = Tyf^fthA + ^ cos^t costal + e + £ sin^sintyl + £ 2 ]A (0)- 



i[Jl + £ 2 sin ft cos tJl + f 2 -f cos £t sin t Jl + f 2 ] A 1 (0)-i sin t Jl + f 2 (cos f «7 (0)-i sin ftCi(O))}, 



= = {-^1 + esin^cos^ + e - fcos^sintyl + ?]A Q (0) + 
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[yjl + £ 2 cos f i cos t y^l + £ 2 +£ sin f i sin t yjl + £ 2 ] (0) -i sin i ^1 + ^ 2 ( -i sin ^C (0) +cos f tC^ (0) ) } . 

(3-16) 

For instance, for the initial state \ip) = ^(\co) + l c i)) = TTjOflS 2)|<7, 0) + |g, 0) |g, 2)), we 
have the initial conditions C o (0) = Ci(0) = ^ and A o (0) = Ai(0) = 0. Then, by using the 

relations (3-16), the evolved state of the system will take the form 

e -i{l+£)t 

|^(t)> = ^^^{[ v / l + ecos^cost v /l + e-^sint v /l + e](|co) + |ci))- ? sint v /l + e(|a ) + |a 1 ))}. 

Now, in order to investigate generation of entanglement between two atoms, we can evaluate 
the density matrix p a associated with the atoms as 

1 



p a (t)=Tr c (W))W)\) 



2(1+ a 



{2(e + cos 2 tJl + e)(\g)W(g\ ® \g) (2) (g\) + 



sin 2 ty/T+?(\e)W(e\ ® \g) {2) (g\ + \g) W (g\ ® |e>t 2 ><e| + |e>< 1 >< £ r| ® |s> (2) (e| + \g)^(e\ ® |e>( 2 )<<?|)} 

where, Tr c denotes the partial trace over the photonic states |2,0), |0,2) and |0,0). Now, for 
a given hopping parameter £, one can use the Peres-Horodecki criteria [24, 25] known also 
as positive partial transpose (PPT) criteria, in order to determine that for which times t, the 
state p a (t) is entangled, and particularly we can obtain the time T at which the perfect transfer 
of photonic entanglement to the atomic one, is achieved. To this end, we choose the order of 
atomic basis as \g,g), \g,e), \e,g) and |e,e), so that the partial transpose of the atomic state 
takes the following matrix form 

1 2(£ 2 + cos 2 t v / T+~F) sin 2 tv/TTF^ 

sin 2 tv/TTF 

sin^v/TTF 

^ sin 2 ty/1 + £ 2 j 

The corresponding eigenvalues of (p a (t)) Tl are given by A = sin 2 ty/1 + £ 2 with double degen- 
eracy, and A± = (£ 2 + cos 2 ty/1 + £ 2 ) ± \J (£ 2 + cos 2 ty/1 + £ 2 ) 2 + sin4 ty/T+^. Therefore, the 



(Pa(t)) 



Ti 



eigenvalue A_ is clearly negative, except for times t 



k-K 



=, k e Z, where the atomic state 
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p a is separable. In order to evaluate the amount of entanglement of the atomic state p a , one 
can calculate the corresponding concurrence [26], as C(p a (t)) = sin2 Vp^ Then> for times 
T = (2 ' +1) " , I G Z, (or in the suitable units, T = (2 ' +1)7r with f = -J^) the maximum 
value of the atomic entanglement is achieved and the corresponding concurrence takes its max- 
imum value C max = Y+^2> where the atomic density matrix will be maximally entangled for 
small hopping £ — > (see Fig.l for different values of £). Moreover, this result indicates that, 
for large hopping strength £ — > oo, we have C(p a (t)) — > and so, the initially unentangled 
atoms, remain effectively unentangled for all the next times. 

Now, consider the initial state |^(0)) = |e,0)|#,0) with C o (0) = Ci(0) = A^O) = and 
A o (0) = 1. The equations (3-16) give 

— ie~^ i 

Cq ^) = 1 1 i c .) SintJl + £ 2 cos£t, 

= li-^ sin^/l + ^sinft, 
V 1 + ? 

—it 

A,(t) = -^= {^1 + e cos ty/l + e cos ft + £ sin t^l + e sin ft}, 

Ax(t) = V 1 + £ 2 cos t^l + ft sin ft - £ sin t ^1 + ft cos ft}, 

so that, one obtains 



law + = ****** 



2 



I . ,, N|2 , . /,m 2 £ 2 + COS 2 VI + £ 2 

\A (t)\ 2 + |^i(t)| 2 = * 

Therefore, in the limit of large hopping ft we have |A (t)| 2 + | (^) | 2 — > 1 and so the only 
effective terms which survive, are A (t) and Ai(t). 

In addition, one can discuss two photon transfer from the first cavity to the second one. 
To do so, we consider the initial state |?/>(0)) = \g,2)\g,0) with initial conditions C o (0) = 1, 
Ci(0) = Aq(0) = Ai(0) = 0. Then, the equations (3-16) give 

C (t) = V 1 + £ 2 cos ft cos tyfi + e - e sin ft sin ty/l + ft}, 
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• —it 

d(0 = ^ { yfl + e sin ft cos t V / T+e + e cos ft sin t ^1 + £ 2 } , 



— ^ / 

= y rT ^ cos^sint v /l + e 2 , 



= ^-^ sinftsint^l + e 



Now, for large enough ^ >, we obtain 

C (*) = e- it cos2^, C^t) = -te~ u sm2&, A (t) = Zll^R^l , A l (t) = ~ & 1^L& . 

Then, after times T = (2fc + lK with k e Z, we have |^(T)) = (-l) fe+1 i e -%, 0)|#, 2) with 
|Ci(T)| 2 = 1 and so, two photons of the first cavity are transmitted to the other cavity, 
perfectly. 



4 Three coupled cavities: Large and Small hoppings 

The case n = 3 identical cavities can be considered similar to the case of n — 2, by using the 
relations (2-7)- (2-9). Here we consider only the limits of large and small hopping ft 

4.1 Large hopping £ ^> 

In the limit of large ft we use the equation (2-12) to obtain 

Co(t) « i{[C (0) + d(0) + C 2 (0)]e- 4i ^ + [2C (0) - d(0) - C 2 (0)]e 2 ^}, 

Ci(t) ~ i{[C (0) + d(0) + C 2 (0)]e- 4 ^ + [-Co(O) + 2d(0) - C 2 (0)]e 2 ^}, 

C 2 (t) « i{[C (0) + d(0) + C 2 (0)]e- 4 ^ + [-Co(O) - d(0) + 2C 2 (0)]e 2 ^}; 

A,(t)«A,(0), for Z = 0,1,2. 

The above results are in correspondence with those of Ref.[18]. By considering the initial state 
|^(0)) = \c ) = \g, 2) \g, 0)\g,0) with initial conditions C o (0) = 1, Ci(0) = C 2 (0) = A (0) = 
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Al(0) = A 2 (0) = 0, we obtain the evolved state after time t as 

\m) ~ ^{(e" 4 * + 2e 2 ^)|c ) + (e- 4 ^ - e 2 ^)(| Cl > + |c 2 ))}. 

Then, the probability of observing two photons at the first cavity (|C (t)| 2 ) and that of ob- 
serving two photons at the two other cavities (\C\(t)\ 2 = |C 2 (t)| 2 ), are given respectively by 

|Co(i)| 2 = ^[1 + 4(1 + cos6£t)], \Ci{t)\ 2 = \C 2 {t)\ 2 = ^[1 - cos6#], 

which indicates that for times T = ( 2fc + 1 ) 7r ; with k E Z, the corresponding two photons initially 
located at the first cavity, are transmitted to one of the other cavities with equal probability 

4 
9' 



Cb(t) 
Ci(t) 

C2(t) 

Mt) 
Mt) 



4.2 Small hopping £ <C 

In the limit of small hopping £ <C, the equation (2-15) leads to the following results for three 
identical cavities: 

{(e~ 3 ^+2) 

{(e" 3 ^+2) 
{(e- 3 ^+2) 
{(e" 3 ^+2) 
{(e- 3 ^+2) 
{( e - 3 ^+2) 

Now, by considering for example the initial state | -0(0) ) = \ao) = \e,0)\g,0)\g,0) , with 
initial conditions A {0) = 1 and C o (0) = Ci(0) = C 2 (0) = Ai(0) = A 2 (0) = 0, we obtain 
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Mt) ~ 3 (e -3 ^ + 2), = A 2 (t) « (e" 3 *' - 1). 

Therefore, after times t = /c7r, with k E Z, the probability amplitudes C (t), Ci(i) and C 2 (t) 
will be zero and we will have A (kn) « (-i)*^ 1-0 ** ( e -3ifl»r + 2 ) an( i Ai(A;7r) = A 2 (A;7r) « 

6 3 ( — - — ( e - 3l £ fe7r _ i). This indicates that, by choosing the hopping strength as £ = 
with large k & Z and small / G Z, the excitation of the first atom located at the first cavity, 
can be transmitted with equal probability I, to one of the other two atoms. 



5 Conclusion 

In summery, the quantum entanglement properties of n coupled atom-cavity systems via two- 
photon exchange interaction, was analyzed. By employing the photon number operator sym- 
metry of the Hamiltonian, the corresponding Hilbert space's dimension was reduced from 3™ to 
2n and then by introducing some suitable Fourier transformed basis states for the state space of 
the system, the corresponding Hamiltonian was block-diagonalized with 2 dimensional blocks. 
Due to this useful reduction, the corresponding Shrodinger equation was solved exactly for any 
number n of atom-cavities, and excitation and photon transition between the atoms and the 
cavities was discussed. The perfect transfer time and the times at which the maximal quantum 
entanglement can be periodically generated between the initially unentangled atoms are ob- 
tained in terms of the hopping parameter (coupling strength) between the cavities. The large 
and small hopping limits was discussed where, it was shown that for large hopping strength, 
the initially unentangled atoms remain effectively unentangled forever. 
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Figure Caption 

Figure. 1: Shows the concurrence C(p) of the atomic state p a (t) with respect to time, for 
different values of hopping strength (a) £ = 0.1, (b) £ = 0.5, (c) £ = 0.9 and (d) £ = 2. 



